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We have calculated the low-field magnetic susceptibility x °f a system consisting of non- 
interacting mono-dispersed nanoparticles using a classical statistical approach. The model makes 
use of the assumption that the axes of symmetry of all nanoparticles are aligned and oriented at a 
certain angle ip with respect to the external magnetic field. An analytical expression for the tem- 
perature dependence of the susceptibility x(T) above the blocking temperature is obtained. The 
derived expression is a generalization of the Curie law for the case of anisotropic magnetic particles. 
We show that the normalized susceptibility is a universal function of the ratio of the temperature 
over the anisotropy constant for each angle ip. In the case that the easy-axis is perpendicular to the 
magnetic field the susceptibility has a maximum. The temperature of the maximum allows one to 
determine the anisotropy energy. 
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I. INTRODUCTION 

The theoretical investigation of magnetic single- 
domain particles goes back to the pioneering work of 
Stoner and Wohlfarth (lj. In the last decade, the in- 
terest in fine magnetic particles is growing, because of 
the prospect of application in high-density data storage 
devices (see for example Ref. &). The size of the mag- 
netic domains depends on the particular material and 
preparation method. However, particles with dimensions 
less than one micrometer, i.e. nanosized particles, usu- 
ally consist of a single magnetic domain Q. Therefore, 
the term "magnetic nanoparticles" is normally used for 
fine single-domain particles. 

The remarkable properties of magnetic nanoparticles 
are due to their anisotropy. The anisotropy has its ori- 
gin in the non-spherical shape of the particle (shape 
anisotropy) or in spin-orbit interaction. In the most 
simple model, the anisotropic nanoparticle can be rep- 
resented by an ellipsoid of revolution. In this case, the 
anisotropy is described by one single parameter K, which 
is called the anisotropy constant. This constant K can 
be calculated exactly from the ratio of the length of the 
minor and major axis of the ellipsoid 0]. At the block- 
ing temperature Tf, = KV/25kB (V is the volume of the 
particle and ks is Boltzmanns constant) the nanoparti- 
cle shows non-equilibrium properties on the time scale 
of the order of 100 seconds M. In principle, magnetic 
particles are suitable for recording only at temperatures 
below ~ Th/2. This is why a large number of theoretical 
and experimental papers have been devoted to the inves- 
tigation of magnetic nanoparticles at temperatures below 
Ti,. However, up to now only a few theoretical investi- 
gations have been reported focusing on the properties of 
nanoparticles at temperatures above Tf, 



It is well known that the temperature dependence 
of the low-field susceptibility x(T) of a system of 
non-interacting nanoparticles with randomly orientated 
axes of symmetry follows the Curie law ||. In this 
case, x{T) does not contain any information about the 
anisotropy of the nanoparticles. However, the time- 
averaged magnetization of an isolated particle is obvi- 
ously anisotropic. The susceptibility of a system of non- 
interacting nanoparticles with parallel axes of symmetry 
was analytically calculated only for several special cases 
The purpose of the present paper is to report on 
the temperature dependence of the low-field susceptibil- 
ity in the general case of an arbitrary angle between the 
symmetry axis and the external magnetic field. 



II. MODEL 

Let us consider a system consisting of a large number 
of identical non-interacting nanoparticles with magnetic 
moment fi and volume V and all axes of symmetry paral- 
lel to each other. Since the interactions are neglected this 
many-particle problem can be reduced to that of a single 
particle. We assume that the particle has the shape of 
an ellipsoid of revolution with an anisotropy constant K. 
The spin-orbit contribution to the anisotropy is assumed 
to be zero. For a fixed orientation of the particle in space, 
the axis of symmetry makes an angle ip with respect to 
the external magnetic field H. The potential energy of 
the particle is then given by 



E = -KVcos 9 - [iHcosuj, 



(1) 



where 9 is the angle between the magnetic moment and 
the applied field and u> is the angle between the axis of 



1 



symmetry and the applied field (see Fig.l). The first 
term in Eq. (|]) yields the anisotropy energy, while the 
second term represents the Zeeman energy of the mag- 
netic dipole in the external field. Our assumption of non- 
interacting particles is valid if the energy of the dipole- 
dipole interaction /ij/Uj/rf- and \ij are magnetic mo- 
ments of two neighboring particles and is the inter- 
particle distance) is much smaller than all other ener- 
gies occurring in the problem, namely the anisotropy en- 
ergy KV, the Zeeman energy \xli and the thermal energy 
k B T. 
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FIG. 1. Spherical coordinate system used to calculate the 
susceptibility of a single magnetic nanosized particle. 6 is 
the angle between the magnetic moment /i and the axis of 
symmetry (dashed vertical line), u> is the angle between the 
magnetic moment and the field direction H, tp is the angle 
between H and the axis of symmetry, and £ defines the angle 
between the projections of /i and H on the equatorial plane. 



For K positive the energetically favorable direction of 
magnetization (at zero field) is parallel to the symmetry 
axis. Hence, the case K > corresponds to the"easy- 
axis" type of particle, while the case K < represents 
the "easy-plane" type of particle. The difference between 
these two cases is shown in Fig. 2. The prolate shape 
of the particle results in a positive anisotropy constant. 
On the contrary, the oblate shape gives rise to a negative 
value of the anisotropy constant K. The theory pro- 
posed here enables one to describe both cases, however, 
we concentrate mainly on the "easy-axis" type of parti- 
cles. Since "easy-plane" particles do not show blocking of 
the magnetic moment at any temperature, these parti- 
cles cannot be used for recording. The magnetization of 
"easy-plane" nanoparticles will be briefly discussed in a 
separate subsection. 





K>0 



K<0 



FIG. 2. Schematic representation of magnetic nanopar- 
ticles with opposite signs of anisotropy constant. The 
anisotropy is due to the ellipsoidal shape of the particle. 
The prolate shape (a) gives rise to a positive value of the 
anisotropy constant K. The oblate shape (b) gives rise to a 
negative K. tp is the angle between the axis of symmetry and 
the external magnetic field H . 



Above the blocking temperature the magnetic 
nanoparticle is in an equilibrium state, and, therefore, we 
can use a classical statistical approach and classical ther- 
modynamics to calculate the magnetic properties. Thus 
our results are of relevance above Tj,. 

By introducing the parameters a — KV j (fcsT) and 
(3 = fiH/{kBT), the partition function for a single do- 
main particle can be written as |l: 

Z = — / / exp(acos 2 6> + /? cos cu) sin 9d9d£. (2) 
4tt Jo Jo 

The angle u can be expressed through ip and the inte- 
gration variables 8 and £ by the relation: 



cos u> — cos 9 cos ijj + sin 9 sin ip cos £. 



(3) 



Once the analytical expression for the partition func- 
tion Z{a, P, ip) is known one can calculate the component 
of the equilibrium magnetization Mm parallel to the field 



= 



1 dZ 
Mil = a— — . 

11 ' Z df3 



(4) 



In order to obtain the analytical expression for 
Z(a,/3,ip), the double integral in Eq. (^) should be cal- 
culated first. By reducing the double integral for the 
partition function in a single integral M it follows: 



Z 



n/2 



exp(cv cos 2 6*) cosh(/3 cos 9 cos ip) 
x I (pj sin 9 sin tp)sin9d9, 



(5) 
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where Iq(-) is the modified zero-order Bessel function of 
the first kind. The same expression was also derived as 
an intermediate result in Ref. fLOfl . We used Eq. (||) as 
the starting point for our calculations. 



III. THE PARTITION FUNCTION 

By inserting the expansions for Iq(-) and cosh(-) into 
Eq. dE|), and by switching the order of integration and 
summation, we obtain the partition function as a infinite 
series (see Appendix A): 



2d 



n!r(n + 3/2) 

n— 



Q n (a,ip), 



(6) 



where Q n (a,ip) represents the internal sum: 



k=0 



k\{n-k)\ 



cos 2fc ip ■ sin 



2(n-fc) 



(7) 



Here M (■,-,■) is the Kummer function (confluent hyper- 
geometric function) and T(-) is the gamma function. The 
expressions @ and (Q) for the partition function are sim- 
ilar to the ones obtained in Ref. Jl(J , where the non-linear 
magnetisation of superparamagnetic particles with a ran- 
dom orientation of easy axes is analysed. Notice, that 
the function Q n (a,ip) does not contain the parameter (3. 
This significantly simplifies the calculation of the deriva- 
tive dZ/d/3: 



2n-l 



dp £j(n-l)!r(n + 3/2) 



(8) 



Both series in the expressions for Z and dZ/d(3 con- 
verge at all values of a, /3 and ip. The reminder of the 
iV-th partial sum for the partition function and its deriva- 
tive, denoted by Rn+i{Z} and RN+i{dZ/d(3}, respec- 
tively, can be estimated as: 



R 



N+l 



{Z}< 



(0/2) 



2iV 



[(tf+i)!] a 



exp(|a| + -), 



rdZ, (0/2) 2N+1 



I 0/9 J 



(TV!) 



h2 



exp (|c 



(9) 



(10) 



These two inequalities are important for estimating the 
errors in the numerical calculations. 



IV. LOW FIELD SUSCEPTIBILITY 

In the low-field case fiH <§; fc^T ((3 <§; 1), only the 
leading terms in (j^) and (||) are retained. This simplifi- 
cation makes it possible to derive an analytical expression 
for the low- field susceptibility (see Appendix B): 
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FIG. 3. The normalized susceptibility X" ~ X ' KV/ ft 2 as 
a function of the dimensionless temperature t = ksT / (KV) 
for K > at several values of the angle ip between the easy 
axis and the direction of the magnetic field, as indicated. The 
vertical dotted line at tb = 1/25 represents the blocking tem- 
perature, below which the presented model is no longer valid. 
The dashed line corresponds to the Curie law. 
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M(i,l,a) 



(11) 



Notice, that Eq. ([11]) represents the component of the 
magnetic susceptibility parallel to the magnetic field per 
particle. To obtain the volume susceptibility the last 
equation should be multiplied by the density of particles. 
Eq. dll| ) is a generalization of the Curie law for the case 
of an anisotropic magnetic particle. We stress that the 
function to calculate the magnetic susceptibility derived 
here is an analytical expression. Although the Kummer 
function M (•, •, •) is not frequently used, its properties are 
well known. The confluent hypergeometric functions are 
available in most mathematical software packages. The 
validity of Eq. (|ll| ) can easily be checked for two special 
cases: zero anisotropy (K = 0) and randomly orientated 
symmetry axes. In both cases, it gives the Curie law ex- 
pression x = M 2 /(3fcsT). This is in agreement with the 
results presented in Ref. ||. 

In order to compare the calculated susceptibility with 
the result obtained by Chantrell et al. we use the 
integral representation of the Kummer function JTlj ] and 
transform Eq. mm into 



k R T 



{ COS ^ + ^- 3COs2 ^( 1+ 2^-^))}' 



(12) 
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where 1(a) — J_, exp(ax 2 )dx. Our Eq. (|l^ ) differs from 
the result given by Chantrell et al. in the sign of the 
l/(2a) term (see Eq. (7) in Ref. J?)]). Notice, the temper- 
ature dependence of the susceptibility was not examined 
in that paper. 



V. DISCUSSION AND ANALYSIS 



A. Case of positive K 



Let us introduce the new dimensionless quantity 

KV 

Xn = X— o-, 
M 



(13) 



which we call the normalized susceptibility. According to 
Eq. (HI), Xn is a function of the two parameters a and 
tp. The anisotropy constant K and the temperature T 
enter in the expression for \n only as a quotient through 
the parameter a = KV/ (ksT). Therefore, the general 
behavior of the temperature dependence of the normal- 
ized susceptibility is the same for different nanoparticles 
with different anisotropy constants, but with the same 
value of the angle tp. The effect of the constant K is to 
change the temperature scale only. Therefore, we show 
in Fig. 3 the normalized susceptibility Xn a s a function of 
the dimensionless temperature t = \ksT / (KV)\ — \a~ 1 \ 
calculated for the case K > ("easy axis") at several 
values of the angle tp. For comparison we also show the 
susceptibility for particles with randomly oriented direc- 
tions of easy axes. As we mentioned above, the latter 
susceptibility obeys the Curie law. From Fig. 3 it follows 
that the low-field susceptibility for the parallel (tp = 0) 
and perpendicular (tp = tt/2) orientations of the easy 
axis with respect to the external field show very different 
behavior. For -0 = the susceptibility values exceed the 
ones obtained by the Curie law at all temperatures. On 
the contrary, for tp = tt/2 Xn(t) is always smaller than 
obtained by the Curie law. This qualitative difference 
is not very surprising and can be explained in a simple 
picture as follows. If the field direction coincides with 
the easy axis (tp = 0) the anisotropy forces the magnetic 
moment to align along the external field. As a result the 
averaged component of the magnetic moment projected 
on the field direction is larger. In the perpendicular con- 
figuration, i.e. the external field perpendicular to the 
easy axis (tp = tt/2), the magnetic moment of the par- 
ticle is deflected from the direction of the magnetic field 
by the anisotropy, which results in a smaller averaged 
magnetization. 

A quite important result is that Xn(t) goes through a 
maximum at t max = 0.23 in the case ip = tt/2. The value 



of t max corresponds to the temperature: 



Using Eq. (|1J), the anisotropy energy KV can be de- 
termined in a simple way from the experimental curve 
of the temperature dependence of the low-field suscepti- 
bility. However, the maximum becomes weaker and fi- 
nally disappears, when the angle ip deviates from tt/2, 
i.e. when the easy axis is not exactly perpendicular to 
the field. For example, for tp = 97r/20 this deviation 
is only 9°, but Xn(t) no longer exhibits a maximum, as 
shown in Fig. 3. 
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FIG. 4. The normalized reciprocal susceptibility 

X^ 1 = X 1 ^ 2 /(KV) as a function of the dimensionless tem- 
perature t = ksT / (KV) for K > at several values of the 
angle tp between the easy axis and the direction of the mag- 
netic field, as indicated. The dashed line corresponds to the 
Curie law. 



The temperature dependence of the reciprocal dimen- 
sionless susceptibility Xn 1 ^) is shown in Fig. 4, for the 
same values of the angle tp as used in Fig. 3. In the 
high temperature region the reciprocal susceptibility fol- 
lows a Curie- Weiss law x" 1 = 3(t — 9d). The dimen- 
sionless Weiss constant 6d is a function of the angle tp. 
The sign of 0d is positive for tp = and negative for 
tp = tt/2. The Curie- Weiss law and the angular depen- 
dence 9d(tp) can be evaluated from Eq. (|ll|). The case of 
high temperatures corresponds to small values of the pa- 
rameter a. Using an expansion of the Kummer function 
(see Eq. (13.1.2) in Ref. |ll|) the reciprocal susceptibility 

A n 



1 for the case a < 1 can be written as 



X^ = - + ^(l-3cos 2 V<), 
a 5 

from which it immediately follows 



T„ 



0.23KV/k B - 



(14) 



^ = -^(l-3cos 2 ^). 



(15) 



(16) 
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In practice, samples consist of nanoparticles with a sig- 
nificant volume distribution. Moreover, the easy axes are 
usually randomly orientated. As we mentioned above, 
randomly orientated easy axes impede the observation 
of anisotropy at temperatures larger than the blocking 
temperature. A significant volume distribution casts se- 
rious doubts on the applicability of simple models used to 
examine the measured susceptibility. We are quite con- 
fident, however, that the current efforts in improvement 
of the sample preparation process of nanoparticles with a 
narrow size distribution, as well as the development of a 
technique for magnetic orientation of nanoparticles [|l2| , 
will enable experimental tests of the model proposed here 
in the near future. 
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FIG. 5. The normalized reciprocal susceptibility 

X^ 1 = X~V 2 /{KV) as a function of the dimensionless tem- 
perature t — }zbT / (KV) for K < at several values of the 
angle i[>, as indicated. The dashed line corresponds to the 
Curie law. 
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B. Case of negative K 



APPENDIX A: 



The model proposed here can also be applied to cal- 
culate the low-field susceptibility of "easy-plane" type 
particles. In Fig. 5 we show the calculated temperature 
variation of \n 1 W f° r several values of ip. Unlike the case 
of K > 0, the magnetic susceptibility of an "easy-plane" - 
type particle is a monotonic function of temperature for 
all tp. Such a particle does not exhibit blocking of the 
magnetic moment, because states with minimal poten- 
tial energy are not separated by an anisotropy barrier at 
zero field. 



VI. CONCLUDING REMARKS 

In this paper we have presented a simple model to cal- 
culate the susceptibility of nanosized particles above the 
blocking temperature. The model describes the case of 
monodispersed nanoparticles with a perfect alignment of 
the easy axes. We have evaluated the susceptibility for 
arbitrary values of the angle ip, i.e. the angle between the 
field direction and the easy axis. One of the main results 
of our calculations is that the normalized susceptibility 
Xn versus the dimensionless temperature t is a universal 
function for each ip, which does not depend on the par- 
ticular value of the anisotropy constant K. One of the 
advantages of the model is that only one phenomenolog- 
ical parameter is involved, namely the anisotropy energy 
KV. 



Using expansions for Iq(-) and cosh(-) (see Eq. (9.6.12) 
and Eq. (4.5.63) from Ref. Jll|]) the partition function 
can be written as: 



Z(a,f3,ip) 



tt/2 



exp(a cos 0) > - 



2i 



i=0 



(2i) 



^ (I/3sin(9sinV>) 2j 
x > ^ttt; sin t 

3=0 



if- 



;H2 



(Al) 



By changing the order of summation and integration and 
using relation (2i)l = 2 2z iW(i + 1/2)/y / 7t , which can be 
obtained from Eq. (6.1.12) in Ref. [ll| we arrive at: 

Z(a,/^) = V^ — -— - ^ 



^2«iir(i + i)^ (j!) 



H2 



,7T/2 

/ exp(acos 2 60cos 2l 6>sin 2j+1 6>d6>. (A2) 
Jo 



The integral in Eq. (A2) can be transformed to the inte- 
gral representation of Kummer function (see Eq. (13.2.1) 
in Ref. |ll[) by substituting t = cos 2 8. Then the parti- 
tion function is: 

^ (l/3cos^) 2 ^(l/3sin^ 

i=0 ' j=Q J ' 

x^±i^±M. (A3) 

r(i + j + §) 
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By changing the order of summation (k = i, n = i + j) 
to gather terms with the same power of the parameter f3 
we obtain: 



By substituting the expressions for f2o and Qi, ob- 
tained from Eq. (]?]), and taking into account that 
Af (£,_£, a) = 3M(i |,a)-2M(i f,a) (see (13.4.3) in 



(i/3cosV) 2fc (i/3sin^) 2 ("- fc ) 

n=0 fc=0 



Ref. 



' 2 1 



we arrive at the following expression for My 



fc! 



(n-A;)! 



M(fe+i n+|,a) 



r(n + 4) 



(A4) 



Mi, = npfcos 2 il>+(\- cos 2 M<y2 ' 2 ' 



2 ' 2 ' 



(B4) 



Last equation can also be written as: 



E 



^nir(»+|)j^*I(»-*)! 

xcos 2fe i/' sin 2( "- fe) r/> • M(k + -, n + |, a). (A5) 

Since multiplying the partition function by any constant 
does not affect th e m agnetization the prefactor y/n can 
be omitted. Eq. QAq) is identical to the expression for 
Z(a,f3,ijj) in Eq. (|6|), except for the prefactor yfn. 



APPENDIX B: 

Considering only the leading terms in (^|) and (||) we 
write: 



Z(a,P,<ip) 



r(3/2) ' 



dZ _ 13 fii(a,ip) 
d/3 ~ 2 ' r(5/2) ' 



(Bl) 



(B2) 



The parallel component of the magnetization in low 
magnetic field is then given by: 



Mil = U ; -. 

11 ^3fi (a» 



(B3) 



From this last expression for Mm we obtain Eq. ( |11| ) using 
the relationships j3 = afiH / (KV) and \ — Mii/H. 
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